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Abstract: It is shown that using elliptically polarized light permits selecting 
well-defined subsurface volumes in a turbid medium. This suggests the 
possibility of probing biological tissues at specific depths. First, we present 
the method and preliminary results obtained on an Intralipid phantom. We 
next report on the method's performance on a biological phantom (chicken 
breast) and, finally, on the exposed cortex of an anesthetized rat. 
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1. Introduction 

2D optical imaging has proven a valuable tool for functional exploration in biological research 
such as in neuroscience (imaging the exposed cortex), ophthalmology (retinal imaging), or 
oncology (skin, uterus cancers detection and monitoring) [1-4]. Yet, there remains concern 
over the fact that images are two-dimensional (2D) projections of a combination of signals 
from superficial and deeper tissue layers, with an a priori unknown proportion. Conventional 
optical imaging of intrinsic signals (IOS) provides high-resolution functional 2D images, by 
monitoring changes in absorption of intrinsic chromophores. However, such systems do not 
allow depth discrimination. As in vivo IOS is usually performed under reflectance geometry 
[5], illumination sources and detectors belong to the same half-space. Under these conditions, 
these imaging systems suffer from (i) blurring due to mirror reflections at the interface; (ii) 
loss of depth information about the signals' origin, due to an indiscriminate collection of all 
photons, irrespective of their pathlengths. To tackle the first problem, one can use linear 
polarization illumination with cross-polarization detection, but in that case, the unpolarized 
backscattered detected photons are those who have undergone a large number of scattering 
events and hence a large proportion is likely to come from deep locations. Moreover, the 
selective detection of photons stemming from the surface or other specific depths remains a 
challenge. 

The technique exposed hereafter is based on a polarimetric selection of the photons. 
Intuitively, when travelling through the medium, polarized photons will maintain their state of 
polarization to a different extent, depending on the absorption or scattering events, that is, 
according to the transport mean free path (MFP) [5-9]. The use of polarizing filters allows, up 
to a point, selecting photons for which polarization state is maintained: this is the basis of the 
polarization gating imaging methods. Linear polarization filters have been first extensively 
employed in sub-surface imaging, where one aims for surface mirror reflection elimination. A 
widely used technique is then detection via the cross-linear imaging channel to get rid of the 
mirror reflections. For tissues examination, this has initially been suggested by Anderson [10]. 
Applied to skin examination, these measurements are sensitive to depths superior to 300 |xm 
[11]. Surface contrast enhancement [12,13] is the other main polarization gating application. 
A simple subtraction between collinear and cross-linear imaging channels allows the 
separation of the surface (or shallow subsurface) contribution from the multiply scattered deep 
volume part, resulting in an enhancement of the surface image. This has also been applied to 
the determination of the optical properties of biological tissues [14]. Polarization gating is also 
used for localized spectroscopy [15,16] in order to extract surface signatures. Ellipsometry has 
also been performed on biological tissues: linear polarized illumination has been coupled with 
tunable polarization detection in order to characterize anisotropic structures in tissues, such as 
collagen fibers [17]. 

However, linear polarization gating, alone, does not allow filtering at the same time mirror 
reflections and multiple scattered light. In that sense, the technique proposed by Morgan and 
Stockford [18], advocating the use of circularly polarized light, is very appealing as it offers 
the possibility of filtering both mirror reflections and multiple scattered light, without any 
special preparation of the samples [11,19]. Since biological tissues are Mie scatterers [5], this 
filtration with circularly polarized light can indeed be achieved in tissues, allowing a deeper 
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screening than linearly polarized light. MacKintosh et al. [6] introduced the concept of 
"polarization memory", describing the way the initial polarization is maintained through the 
variation of the number of scattering events, and showed that the polarization of circularly 
polarized light is indeed maintained through a larger number of scattering events than that of 
linearly polarized light. This effect has been demonstrated in biological tissues and used for 
the localization or analysis of specific sample regions not only in transmission image 
acquisition geometry [20] but also in reflection geometry [21], which is more adapted to 
biological tissue screening in vivo. 

Here, we adopted the approach proposed by [18]. In addition, we introduce the novel idea 
that it is possible to select the probed depth as a function of the polarization ellipticity, thus 
allowing screening the tissue at any specific depth comprised between the surface and a 
maximum depth defined by the maximum penetration depth of the circularly polarized 
photons. An experimental setup has been developed and the principle is demonstrated through 
phantom experiments. Our final motivation in the present work is to test the feasibility of the 
technique for depth resolved IOI for exposed cortex in vivo examination. A first experiment 
has been conducted on an anesthetized rat and is reported as a first illustration. 

2. Materials and methods 

2. 1. Polarization gating for volume selection 

The principle has been developed by Morgan et al. [18]. It consists in sequentially 
illuminating the biological tissue with light polarized linearly or circularly. Detection 
performed under cross-linear polarization allows the collection of multiple scattered photons 
(those coming prevalently from deep tissues), whereas co-linear detection is restricted to the 
collection of linear polarization maintaining photons (mostly surface probing). If circularly 
polarized light is considered, one can distinguish three state of polarization domains: (i) the 
first backscattered photons have undergone a mirror reflection and thus flipped their helicity; 
when analyzed, those contribute to the cross-circular component (thus usable for the 
validation of the surface, mirror reflected, image); (ii) the second kind of backscattered 
photons emerge after a series of forward scattering events and maintain their original 
polarization state. They contribute to the co-circular component; (iii) the third type of photons 
have been randomly polarized by multiple scattering events and contribute equally to all 
polarization states. 

By subtracting cross-linear images — to which the contribution of multiple scattered 
photons is predominant — from co-circular polarization images, one can probe a medium- 
depth subsurface volume. The latter is free of mirror reflections, while eliminating the 
multiple scattered photons coming from greater depths (probed by the cross-linear image). 
The extent of this middle-depth, "subsurface" volume depends on the optical properties of the 
medium. 

Importantly, the same reasoning as above can be applied to elliptical polarization. 
However, the probed volume will be shallower than for circular polarization. 

2.2. Experimental setup and image processing 

The experimental setup (Fig. 1) is composed of a 150W halogen light lamp (OSL1-EC, 
Thorlabs SAS, France) for illumination. The beam is homogenized through a telecentric 
objective, and wavelength filtered with an interference chromatic filter (632.8 nm, lOnm 
FWHW, Thorlabs SAS, France). The detector is a CCD camera (12 bits, Orca ER, 
Hamamatsu Photonics, Japan) with an objective (7000E, Laser Components SAS, France) for 
detection. The incident beam passes through a polarizer and quarter wave plate. In the 
detection path there are also a quarter wave plate and an analyzer. The relative rotation angle 
between polarizers and quarter wave plates is tuned such as to obtain the desired polarization 
modes for illumination and detection. 

The axes of the different optical components were first set with calibrated measurements. 
Linear polarization is obtained, for a given polarizer/quarter-wave plate couple, when the axis 
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Fig. 1. Sketch (left) and photograph (right) of the experimental setup. 



of the polarizer and the quarter- wave plate are aligned; circularly (right or left) polarized light 
is obtained when the fast axis of the quarter- wave plate is 9 = 45° (+ or -, respectively) off the 
linear polarizer axis. 

Experiments have been performed under different illumination/detection channels: (i) 
"collinear" image (illumination under linear polarization and detection under co-linear 
polarization) detects surface mirror reflected signal plus volume scattered signal; (ii) 
"crosslinear" image (illumination under linear polarization and detection under cross-linear 
polarization) detects only volume scattered signal; (iii) "co-elliptical" images (illumination 
under elliptical polarization and detection under co-elliptical polarization) detect subsurface 
signal and volume scattered signal. 

With these schemes, the following imaging channels can be defined as follows: 

(A) "surface" = "collinear" - "crosslinear"; 

(B) "middle" = "co-elliptical" - "crosslinear"; 

(C) "multiple scattered" = "crosslinear". 

Evaluation of the polarization-maintained signal from the middle channels is based on the 
assumption that the fraction and characteristics of the multiple scattered photons are the same 
for linearly and elliptically polarized light. This approximation is substantiated in [19] for the 
case of circularly polarized light. In this reference, a Monte Carlo simulation shows that the 
penetration depth profiles of the collected multiple scattered photons are the same under 
cross-circular and cross-linear configurations, even if the number of multiple scattered 
photons collected under cross-circular configuration is smaller than the number collected 
under cross-linear configuration, which is coherent with the above-mentioned "polarization 
memory" effect. The image resulting from the subtraction of the middle channel and the 
multiple scattered channel, approximated by the measurement under cross-linear 
configuration, contains indeed only polarization-maintained signal, but this signal is 
underestimated, leading to images with a poorer signal to noise ratio (SNR). Although the 
reasons of the similarity of the depth distributions of multiple scattered photons, collected 
under both polarization configurations, are not obvious, the same approximation has been 
adopted here for light polarized elliptically. A numerical verification (Monte Carlo 
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simulations not shown here) has shown indeed this approximation has no major impact on the 
resulting images, provided, at least, the anisotropy factor is high, which is the case in 
biological tissues. The validity of this approximation depends indeed on the optical properties 
of the medium and deserves a specific study, which is out of the focus of this paper. 

The exposure time of the camera was adapted to each measurement in order to preserve a 
high SNR for each measurement. Yet, to facilitate comparison between different images, the 
data are expressed per unit exposure time. No other processing has been applied to the images. 

2.3. Samples description 

Experiments have been conducted on three kinds of samples: the first series on Intralipid 
liquid phantoms for quantification evaluation, the second on a piece of chicken breast for 
biological tissue feasibility, and the third one, on an anesthetized rat with exposed cortex, as 
in vivo test. 

2.3.1. Liquid phantom 

The liquid phantom was composed of aqueous Intralipid (20%, Sigma-Aldrich, France) 
dilution at adequate concentrations. A target in the form of a ruler (plastic) had been placed 
obliquely in a tank according to the sketch in Fig. 1 . Placed in that way, each graduation read 
on the CCD screen corresponded to a different depth. The concentration of Intralipid had been 
adapted to match biological tissues optical properties but with a scaling in the size of the 
sample for macroscopic measurement of depth. The optical properties of Intralipid 1% have 
been determined according to the integral reflectance method [22] : the absorption coefficient 
is negligible and a n' s i% = 10.3 ± 0.5cm -1 reduced scattering coefficient value was estimated. 
According to the accuracy of our optical (magnification) and mechanical components, a 0.1% 
diluted solution had been adopted, corresponding to fi' s o.i% = 0.95 ± 0.5cm" 1 , allowing depth 
measurements at a millimetric scale. The anisotropy factor of Intralipid has been estimated to 
g = 0.73 [23] leading to a scattering mean free path MFP = {\-g)lfi' s ~2.% mm (that is 
approximately 10 times longer than in biological tissues). 

2.3.2. Biological phantom 

The second phantom was a piece of chicken breast. A graphite stick (diameter, ca. 1mm) was 
inserted with a small tilt angle (exact value not evaluable) to serve as inclusion. Caution was 
taken though to keep the superficial tissue sufficiently hydrated. This helped tracking the 
mirror reflections. 

2.3.3. In vivo experiment 

After initial anaesthesia with 5% isofluorane, an adult male Wistar-Han rat was anesthetized 
with Urethane (intraperithoneal injection, 130 mg/lOOg body weight) placed in a stereotaxic 
frame for mechanical stability during craniotomy and during all the duration of the recordings. 
Body temperature was monitored, and maintained between 36 and 38°C with the use of a 
warm water bag. A craniotomy was performed creating a cranial window of ca. 5mmx5mm 
above the barrel cortex. To prevent the cortex from drying, a chamber of dental cement was 
created around the entire cranial opening such as to create a chamber that was filled with 
agarose in saline (2%). To provide a flat optical interface, the chamber was sealed with a 
microscope coverslip. In compliance with the directions of the Local Ethics Committee, we 
kept the number of animals used for this study as small as possible. All procedures were in 
agreement with NIH Guidelines. 

3. Results 

The two phantoms described above have been examined under the various described imaging 
channels. Hereafter, we show the resulting images obtained for the different "middle" 
channels and the "volume" channel. 
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The results obtained in Intralipid 0.1% are represented on Fig. 2. In the top part of the 
image, the different images are displayed with the same colorbar and show that the ruler can 
be seen in progressively larger depth, with a signal whose amplitude increases as well, 
showing that an increasing number of photons are collected as the polarization ellipticity 
increases. The middle part of Fig. 2 displays images with scales adapted (automatic scaling 
between maximum and minimum, with jet color map) to each image in order to enhance the 
contrasts of low signal images (subsurface images). It shows explicitly that some deep 
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Fig. 2. Top: four first rows, "middle" channels images, obtained at different degrees of 
ellipticity, for Intralipid 0.1% experiment. The angle value measures the angle 6 between the 
polarizer and the illumination quarter-wave plate. The same grey level scale was used for all 
images. Fifth row image, deep chamiel image obtained under crosslinear detection, with the 
same field of view but different grey level scale (see colorbar). X axis indicates the actual 
depth of the ruler in the Intralipid, increasing from left to right, while Y axis indicates the 
actual dimension read at the surface of the probed medium. Middle: selected images with color 
level scale adjusted to each image, such as to go from its minimum to its maximum. For all 
images, intensity is expressed in arbitrary units. Bottom: ■ maximum depth (expressed in MFP' 
= 1MFP' = 1.05 cm) visible in a given polarization maintained channel as a function of 

the corresponding angle 8 between the polarizer and the illumination quarter-wave plate, the 
fitting line y = 0.035x (linear regression coefficient R 2 = 0.89) is represented (black line). 
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structures are not visible on subsurface channels while clearly visible on deeper channels: for 
instance, number "7" is not visible on channel 9 = "+8°" while clearly visible on channel 9 = 

"+18°". 

Calculated from the tilt angle of the ruler, here, 1 mm read on the image corresponds to 
0.375 mm depth shift: the horizontal axis on the multiple scattered image represents the 
corresponding actual depth of the ruler in the Intralipid, while the vertical axis indicates the 
actual dimension read at the surface of the probed medium. This experiment shows that 
polarized light can be collected up to ~5 mm in depth (corresponding to 0.5 mm in biological 
tissues) without any difficulty. To fix the ideas, the depth of the last visible bar of the ruler 
(Maximum visible depth), read in a given polarization maintained channel, is plotted as a 
function of the corresponding angle 9. In order to have an estimation of the depth accessible in 
tissues, this depth is reported with a MFP' (transport mean free path) scale, accounting for all 
optical properties. A linear behavior is found showing a complete correlation between the 
polarization ellipticity and the probed depth, up to 9 — "+30°". This experiment can be easily 
reproduced for different optical properties media (for different absorption, scattering 
coefficients and anisotropy factors) and could serve as calibration measurement for tissue 
studies. 

As the result of images subtraction, the SNR of the resulting channel images is smaller 
than that of the measured images. Using basic image processing should allow to reduce the 
noise contain of the images. Interestingly, a specific statistical analysis of the images as the 
Pearson correlation analysis proposed in [13] could improve the accuracy of the estimation of 
the probed depth. 

The chicken breast experiment (Fig. 3) illustrates the feasibility in other kind of tissues. 
The graphite stick had been inserted from left to right with a small angle such that its depth is 
higher on the right hand side of the image. The horizontal and vertical axes, reported only on 
the multiple scattered images, indicate the field of view at the surface of the probed medium, 
in pixels unit, 1 pixel = 39 urn x 39 urn. The depth of the tip of the stick has not been properly 
measured: the coefficient of proportionality between the horizontal axis and the actual depth is 
not reported. However, we reported on the bottom part of Fig. 3 the maximal position (pixel 
units) of the visible part of the graphite stick, read on the horizontal axis of images in Fig. 3, 
as a function of the angle 9. Basically the same comments as above could be done, small angle 
9 representing near surface images, that is small polarization ellipticities, show only the 
superficial part of the stick (left side of the image) while for larger angle 9, that is higher 
polarization ellipticities, right side becomes more and more visible, with an increase of the 
signal intensity. 

The in vivo results are reported on Fig. 4. A network of blood vessels having different 
sizes can be seen on the white light image (Fig. 4, Top left). We focused on a portion of the 
image (Fig. 4, Top right) showing two distinct absorbing structures of different sizes: a very 
large vessel (a), and a smaller one (b). The first one (a) is a large superficial structure and is 
present in all channels. The smaller vessel (b) is visible only beyond a certain degree of 
ellipticity (9 > +14°). Moreover, the signal intensity increases more above vessel (a) than 
below it. This discrepancy is not visible in the deep channel image, suggesting that the 
observed local absorption pattern is restricted to a more superficial volume. 

4. Conclusion 

Here, we reported on the usage of polarization gating imaging with elliptically polarized light. 
Experiments were conducted both on phantoms and in vivo. Our results show that using 
imaging channels with elliptically polarized light allows probing subsurface volumes, which 
are well defined. Whereas the exact depth and extent of these subsurface volumes depend on 
the optical properties of the medium, such polarization gating allows producing images free of 
surface and deep volume scattering contaminations. Our phantom experiments show that the 
method can be used to probe tissues in depth up to at least 0.5 mm, which offers the 
possibility of screening a variety of layered biological tissues such as the exposed cortex, as 
tested here, or the skin. Further work, including modelling of light propagation, improvement 
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of the setup and application of signal processing, should allow using the method to image 
selectively sub-volumes at well-defined, specific depths. In this sense, the presented approach 
constitutes a first step towards functional optical imaging at user-selected depths within the 
tissue. 
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Fig. 3. Same as in Fig. 2 for chicken breast experiment. The depth of a graphite stick (arrow) in 
the tissue increases from left to right. Horizontal and vertical axes of images represent the field 
of view. Pixel size: 39 um. 
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Fig. 4. Top, white light image of the imaging window with a large field of view (left), white 
rectangle marks the reduced field of view in which the image processing is performed (right). 
Middle and bottom, same as in Fig. 2, for rat cortex examination. Pixel size: 42 fim. 
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